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We study the particle motion in the space-time of a Kehagias-Sfetsos (KS) black hole. This is 
a static spherically symmetric solution of a Hofava-Lifshitz gravity model that reduces to General 
Relativity in the IR limit and deviates slightly from detailed balance. Taking the viewpoint that 
the model is essentially a (3-f l)-dimensional modification of General Relativity we use the geodesic 
equation to determine the motion of massive and massless particles. We solve the geodesic equation 
exactly by using numerical techniques. We find that neither massless nor massive particles with 
non- vanishing angular momentum can reach the singularity at r = 0. Next to bound and escape 
orbits that are also present in the Schwarzschild space-time we find that new types of orbits exist: 
manyworld bound orbits as well as two-world escape orbits. We also discuss observables such as the 
perihelion shift and the light deflection. 
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I. INTRODUCTION 



Motivated by the study of quantum critical phase transitions Hofava introduced a (3+l)-diniensional quantum 
gravity model, later on called Hofava-Lifshitz (HL) gravity, that is power-counting renormalizable [H, (see also Q 
for a recent status report). This model reduces to General Relativity (GR) in the infrared (IR) limit, i.e. at large 
distances, however breaks Lorentz symmetry in the ultraviolet (UV), i.e. at short distances. The reason for this is 
that the model contains an anisotropic scaling with dynamical critical exponent z of the form 

r-^ bf , t~>bH. (I.l) 

In the IR the exponent becomes z = 1 and the theory is lorentz-invariant. However, in the UV there is a strong 
asymmetry between space and time with z > 1. In (3+1) dimensions z = 3 [2| and the gravity theory becomes power- 
counting renormalizable. Concretely, this model breaks Lorentz invariance at short distances because it contains only 
higher order spatial derivatives in the action, while higher order temporal derivatives (which would lead to ghost 
degrees of freedom) do not appear. 

A number of explicit solutions of HL gravity have been found, in particular spherically symmetric black hole solutions 
[3-0|- T"^^ most general spherically symmetric solution has been given in and rotating generalizations have been 
studied in Q. One of the open problems of the model is how to couple it to matter fields. The question of how to 
describe particle motion in HL gravity, i.e. to find the equivalent to the geodesic equation of GR has been addressed in 
[lol-[l2l| . In [l^ particles were studied as the optical limi^ of a scalar fieM, while in P a super Hamihonian formahsm 
with modified dispersion relations was used. In both papers it was found that new features arise in HL gravity 
such as superluminal motion and luminal motion of massive particles. In [l2j | a particle action preserving foliation 
diffeomorphisms was introduced and it was found that massless particles follow GR geodesies, while the trajectories 
of massive particles depend on their mass. In most studies of test particle motion the hypothetical corrections to the 
GR geodesies were neglected [l3l - l23 |. 

In this paper we take the latter viewpoint and study solutions to the GR geodesic equation in HL black hole 
space-times, in particular in the space-time of a Kehagias-Sfetsos (KS) solution, a static and spherically symmetric 
solution to HL gravity with vanishing cosmological constant. The geodesic motion in this space-time has been studied 
previously |13l424| and a number of constraints on the parameters of HL gravity have been found. Observables such 
as the perihelion shift and the light deflection were also studied in these papers, however, either approximations were 
used or only circular orbits were studied. In this paper we are aiming at solving the geodesic equation exactly by 
using numerical techniques and at exploring the complete set of solutions of the geodesic equation. 

Our paper is organized as follows: in Section II we give the model and the black hole solutions. In Section HI we 
give the geodesic equation, while Section IV contains our results. We conclude in Section V. 



II. THE MODEL 



A. The action 



The model proposed by Hofava [l|, uses the ADM decomposition of the metric that reads as follows 

ds^ = -N'^df + g,j {dx' + N'dt) {dx^ + N^dt) , (II.l) 

where N{t^x'^) and A^*(i,x') are the lapse and shift functions, respectively, and gij{t,x'^) is the 3-metric with i,j = 
1,2,3. In [4] it was assumed that the theory is invariant under space-independent time reparametrization and time- 
dependent spatial diffeomorphisms, i.e. under 



t i{t) , x' x\t,x') , 
which restricts the lapse function to depend only on t. The action proposed in [2!| then reads 

5 = 5*0 + 5*0 + S'l , 

where 
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So = / dtd^x^N 



, 5o = / dtd^x^N 



and 



Si = / dtd^x^N 



32(1 -3A) 2u;4 *J 2w^ J A- 8 ^ 



(H.2) 
(H.3) 
(II.4) 

(H.5) 
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where 

(7 is the determinant of the metric gij and Rij , Vi correspond to the spatial components of the covariant derivative 
and the Ricci tensor, respectively. C*^ is the Cotton tensor and A, k, /x, w and Kw are constants. The integrand of 
— (5*0 + 51) is interpreted as the potential part, while the integrand of Sq is interpreted as the kinetic part. 
In the IR limit the action is dominated by Sq + Sq and reduces to the Einstein-Hilbert action for 



where c is the speed of light, Gm is Newton's constant and A is the cosmological constant. Note that for A > 1/3, i.e. 
in particular for A = 1, the constant A^y and hence A should be negative. In the following we will set A 1 (unless 
otherwise stated) and consider the additional terms of Si as a (3+l)-dimensional modification of General Relativity. 
The action as given above satisfies the requirement of detailed balance which essentially means that the potential V 
in the Hofava-Lifshitz action derives from a superpotential W: 

I, 1 5W 

and C/'J'^' = \ {g^^9''^ + g^^'d''^) — ^9^'' 9^^ is the DeWitt metric. The requirement of detailed balance drastically reduces 
the number of invariants to consider in the potential V . 

The problem with the theory as stated above is that for A ~ 1 it predicts the wrong sign of the 4-dimensional 
cosmological constant. Moreover the detailed balance condition is chosen solely to simplify the theory. Hence, 
theories that violate detailed balance have been considered. In [3] the following term was added to the action 

S^^ j dtd'x^N^^^^^ujR , (II.9) 

where to is an arbitrary constant. In the Aw = limit which we are mainly interested in here the Einstein-Hilbert 
action is recovered in the IR for 



B. Spherically symmetric solutions 

Kehagias and Sfetsos (KS) found a spherically symmetric, static black hole solution to a Hofava-Lifshitz gravity 
model with action S + Sy for Aw = and A = 1. The Ansatz for the metric is 

ds^ = N^{r)dt^ - f-\r)dr^ - {d9^ + sin^ 9dip^) (11.11) 

and the solution reads 

= f = l + OJr^ - Vuj'^r^ + 4:UJmr , (11.12) 

where w = 16/z^/k^ and m is an integration constant. In (isj constraints on the value of ujm^ were found by comparing 
the perihelion shift in the KS space-time with observations in the solar system. It was found that wto^ > 7.2 • lO^^'^ 
for Mercury, wm^ > 9 • 10~^^ for Mars and ajm^ > 1.7 • 10~^^ for Saturn. Moreover, a similar comparison gave 
ujm? > 8 • 10"^" for the S2 star orbiting the supermassive black hole in our galaxy as well as utrn^ > 1.4 • 10^^^ for 
extrasolar planets [l3| . In [131 constraints from innermost stellar circular orbits (ISCOs) for certain black holes were 
considered and it was found that ui ~ 3.6 • 10~^'*cm~^ (in appropriate units). In [3| the light deflection in the solar 
system was used to constrain the parameter. It was found that Ldm^ > 1.17 • 10""'^^ for Earth, wm^ > 8.28 • 10~^^ 
for Jupiter and wm^ > 8.28 • 10~^^ for the Sun. The IR limit of pi.l2p is given by the Schwarzschild solution 
iV^ = / = 1 - 2m/r. The Kretschmann scalar K = Rf^upaR'^"''" reads 



^ 97-2 ) \ dr ) ^4 + ^4 ' 



(H.13) 
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which for small r behaves like l/r^ . Hence the solution possesses a physical singularity at r = [3] and two horizons 
at 



r± — m ± rn? — — (11.14) 

as long as ujm} > 1/2. Note that the corrections from Hofava-Lifshitz gravity now allow for the existence of up to 
two horizons. The extremal solution has u/rn^ = 1/2 and r_|_ = m. The Hawking temperature of black hole solutions 
is given by Th = k/{2t:), where k is the surface gravity that for static solutions is given by 



K 



.2 



1 



4- 

For the KS solution we find 



'''g''{^^9u){^Jgu) ■ (H.15) 



1 u}{r± - m) 

which in the w — )• c» limit tends to the known Schwarzschild result Th = (87rm)~^. Obviously, the extremal solutions 
with r_L = m have Tjj = 0. For more details about the thermodynamics of black holes in Hofava-Lifshitz gravity see 

e.g. m. 



III. SOLUTIONS TO THE GEODESIC EQUATION IN HORAVA-LIFSHITZ BLACK HOLE 

SPACE-TIMES 

For a general static spherically symmetric solution of the form pi. lip the Lagrangian Cg for a point particle reads 



1 dxf" dx"" _1 1 
2-^'"'^^" 2^" 2 



, dt 



1 /dr 



f \dT 



de\ 
d^J 



I - sin^ I 



dip 
d^ 



(nil) 



where e = for massless particles and e = 1 for massive particles, respectively. 

The constants of motion are the energy E and the angular momentum (direction and absolute value) of the particle. 
We choose = 7r/2 to fix the direction of the angular momentum and have 



Using these constants of motion we get 



dr dr 



(in.2) 



dry ^ 
dr / iV2 



(111.3) 



and 



dr\ / 
d^J ~ LlN^ 



E^-V^sir) 



(in.4) 



where V^ff (r) is the effective potential 



(111.5) 



In the following we will consider the KS black hole solution (jII.12|) . The geodesic equation (|III.4|) then becomes 



1 dr 



r dip 



1 dr 



—T +2 -— P(r)=Q(r), 



dip 



(IH.6) 



where 



P(r) = (cer^ + [e - E' + u^Ll)r^ + L^) 



(in.7) 
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and 

Q{r) = [-2euj{E^ - e)r^ - iujme'^r^ + {-2ujE'^LI + /^ujLle + {E^ - e)^)r^ 

- SujmLler^ + 2LI{ujLI - E^ + - AujmLlr + L^] . (111.8) 

For massive particles (e = 1) the order of the polynomials P{r) and Q{r) is 4 and 6, respectively, while for massless 
particles (e = 0) it is 2 and 4. 



Rewriting pil.6l) we find 



V~Vo = ± I = ■ (111.9) 

-P ± v/P2 + Q 



The motion of test particles in KS black hole space-times has been studied extensively before |13l - ll9l |. however, it 
was not attempted to find the complete set of solutions. This is what we are aiming at here. The integral on the 
right hand of piI.Qp cannot be solved in terms of hyperelliptic functions, at least not to our knowledge. However, an 
analytic treatment seems possible in some limiting cases. This will be reported elsewhere [2^ . In this paper we solve 
the geodesic equation (|III.6p numerically. 



IV. RESULTS 



A. The effective potential 



In order to understand which types of orbits are possible in the KS space-time, we first study the effective potential. 
To make contact with the Schwarzschild case we rewrite (jIII.3[) as follows 

(J)'=f-K=ff(r) , (IV.l) 

where £ = E^ — e and 

Kff(r) = V,s{r) - e = (wr^ - ^/cj^r^ + 4wmr) (^e + + ^ . (IV.2) 

For r ^ {Am/ioy/^ this effective potential becomes Vcs{r ^ (4r7i/aj)^/^) « —2m£/r~2mLl/r'^ + Ll/r'^^ which is just 
the effective potential in the Schwarzschild space-time. 

The first point to note is that while for u: ^ oo the potential at r ^ 1 behaves like T4ff(r ^ 1) ~ —2mLl/r^ (this 
is just the Schwarzschild limit), it behaves like Vcs{r ^ 1) « Ll/r"^ for generic to. Hence there is a positive infinite 
angular momentum barrier for both massive and massless test particles which does not exist in the Schwarzschild limit. 
The first conclusion is hence that test particles with non- vanishing angular momentum cannot reach the singularity at 
r = in the KS space-time. Moreover, for the extremal solution with r = ?■+ = rn we find that dVcff{r) / dr\r=r+ = 
and Vcff(r r+) -e. 

On the other hand, for particles without angular momentum Lz = 0, the effective potential is always negative and 
behaves like Ves{r <C 1) ~ —e\/4:Ujmr for small r, while it is equivalent to the Schwarzschild potential for large r: 
Vcs{r > 1) « -2me/r. 



1. Massive test particles 



In Figs. l(a)|l(c)" we show how the effective potential Vce{r) for a massive test particle (e = 1) changes for different 



values of Lz and lo and m = 1. 

It is obvious that the effective potential at large r doesn't change much when decreasing lj from the Schwarzschild 
limit u — oo. Hence, the types of orbits available for large r are very similar to the Schwarzschild case. This is not 
surprising since Hofava-Lifshitz gravity is a gravity theory that is supposed to modify General Relativity at short 
distances, but has no effects on the long distance physics. In comparison to the Schwarzschild case, the effective 
potential possesses a further minimum at small r. This is represented by the curves in Fig|5] In this latter plot, we 
assume that E"^ — e is a parameter that can have all possible values to show that an additional minimum exists, but 




FIG. 1. The effective potential 14ff(r) for a massive ((a)-(c)) and a massless ((d)-(f)) test particle, respectively, for different 
values of u and L^. 




FIG. 2. The values of — e and l/L^ corresponding to the absolute minimum of the effective potential Vc«{r) at small r for 
different values oi to, m — 1 and e = 1. Note that while here we treat i?^ — e as a parameter that can take arbitrary values, we 
should have — e > —1 when looking for zeros of f — Vc«{r). 



keep in mind that to find the zeros oi £ — V^ff (r) we need to require £ > —1. Note that the value of this minimum is 
negative and always smaller than —1. It increases for decreasing w and becomes equal to —1 in the extremal limit. 
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L = 7.00, E = 0.99, m = 1 and e = 1 

z 




r 



FIG. 3. The two regions of the potential for which bound orbits of massive test particles exist. In region 1, we have manyworld 
bound orbits, while in region 2 there exist bound orbits. Here Lz — 7.0, E'^ — 0.9787, m = 1.0, while u> = 0.51 and uj — 10*, 
respectively. The red dotted-dashed line represents the total energy {E^-e) of the test particle. 



(Si = 5.10, m = 1,e = 1 m = 51000.00, m = 1,e = 1 




(a) oj = 5.1 (b) w = 5.1 ■ 10" 



FIG. 4. The values of — e and 1/L1 corresponding to the maximum (thick upper line) and relative minimum at large r (thin 
lower line) of the effective potential VcR(r) are given for m = 1, e = 1, a; = 5.1 (left) and oj = 5.1 • 10* (right), respectively. In 
the dark shaded region (region I) there exist manyworld bound orbits (MBO) and bound orbits (BO), while in the light shaded 
region (region II) manyworld bound orbits (MBO) as well as escape orbits (EO) exist. In region III there are manyworld bound 
orbits (MBO), while there are two- world escape orbits (TEO) in region IV. 



This is clearly seen in Figs l(a)||l(c' 



These observations lead to the following conclusion for the types of orbits possible which have turning points at 
the minimal radius r = rmin close to r = 0: test particles would move on manyworld bound orbits (MBO) with 
''min < T < r,„ax Or On two-woild escape orbits (TEO) with rmin < r < oo but can never reach r = 0. In comparison 
to bound orbits (BO) and escape orbits (EO), respectively, test particles moving on manyworld or two- world orbits 
cross the two horizons in both directions. That this is always the case for orbits with r^^^ close to r = can be 
seen as follows: since Vcs{'r±) — = ~1 and the turning points are given \yy — e = Vcff (r), the value of rmin is 
always smaller than r_ and the value of Tmax is always larger than r-|_ . In the Schwarzschild space-time manyworld or 
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L = 0.00, m = 1 and e = 1 

z 

0.2 F ' ' ' ^ 




r 



FIG. 5. The effective potential Vcs{r) for radial trajectories (Lz = 0) of massive particles (e = 1) in the space-time of a KS 
black hole with m = 1 and different values of lj. 
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TABLE I. Types of orbits of massive test particles in the KS space-time. The thick lines represent the range of the orbits. The 
turning points are shown by thick dots. The horizons are indicated by double vertical lines. 



two- world orbits are not possible: a particle crossing the horizon would always end at the physical singularity at r 0. 
Note that in the KS space-time we also have bound orbits that are comparable to the bound orbits existing in the 
Schwarzschild space-time. The two regions in which manyworld bound orbits and bound orbits, respectively, exist are 
shown in Figl2]for m = 1 and two different values of w. In region 1 we have manyworld bound orbits, while in region 
2 we have bound orbits. The effective potential varies only little in the region 2 at large r when changing oj from 10^ 
to 0.51, while in region 1 at small r it varies strongly. The above results are summarized in the {E"^ — e)-(l/L^)-plot 
(see Figs. [4(a)p(b)| ). 

The shaded region is bounded by two curves, the one at larger E'^ — e representing the maximum of the potential 
and the other one the local minimum of the potential at large r. The dark shaded region with E'^ — s < (region I) 
corresponds to the values of E^ and for which £ — Vcs{r) has four positive real- valued zeros. Hence, there are two 
different types of orbits: a manyworld bound orbit (MBO) as well as a bound orbit (BO). The light shaded region 
with E^ — e > (region II) corresponds to the values of E'^ and for which £ — V^ff (r) has three positive real- valued 
zeros and hence we have a manyworld bound orbit (MBO) as well as an escape orbit (EO). In the white region with 
E^ — e < (region III) £ — Ves{r) possesses two positive, real- valued zeros such that the corresponding orbit is a 
manyworld bound orbit (MBO). Finally in the white region with E^ — e > (region IV) £ — Kff(r) has one positive, 
real- valued zero and the corresponding orbit is a two- world escape orbit (TEO). These results are also summarized 
in Table U 

Note that the orbits existing in this space-time are very similar to the ones in the Reissner-Nordstrom space-time 
[23, Comparing the case for uj = 5.1 with that for cj = 5.1 • lO'*, we observe that the features of the plot do not 
vary much. This is also true for even smaller values of w. 

Massive test particles with Lz = move on radial geodesies with (p — const.. In this case, the minimum of the 
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m=1, (0 = 0.51, E= 0.99 andL =7.00 



m =1, Q) = 10000.00, E = 0.99 and L = 7.00 





(a) 



Region 1, lu = 0.51 



(b) 



Region 1, lo = 10"^ 



m =1,(0 = 0.51, E= 0.99 and L =7.00 



m =1, (0 = 10000.00, E = 0.99 and L = 7.00 





(c) 



Region 2, ui = 0.51 



(d) 



Region 2, lj = 10'' 



FIG. 6. Examples of manyworld bound orbits (MBO) and bound orbits (BO) of a massive test particle (e = 1) with Lz — 7.00, 
E'^ — 0.9787 in the space-time of a KS black hole with m = 1.00 as well as lj = 0.51 (left) and uj = 10* (right). We show 
manyworld bound orbits (MBO) (region 1, top) and bound orbits (BO) (region 2, bottom), respectively. The red dashed circles 
in the plot represent the horizons of the KS black hole. Note that we are plotting two radial periods during which the particle 
moves from r-min to rmax and back again. 



: 1,M = 0.51,E = 1.80andL =4.00 



1,(0 = 1.53, E = 1.80and L =4.00 





; 1.(0 = 50.00, E = 1.80 andL =4.00 




-6 -4-2 2 



(a) 



uj = 0.51 



(b) 



UJ = 1.53 



(c) 



UJ = 50.00 



FIG. 7. Examples of two-world escape orbits (TEO) of a massive test particle (e — 1) with E — 1.8, = 4 in the space-time of 
a KS black hole with m = 1 and different values of u. The dashed circles correspond to the two horizons of the KS space-time. 



effective potential is at r = tq such tliat 



ro = (£)'^' and V,s{ro) ^ - {2^771^)'^' . (IV.3) 
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positive zeros 


range of r 


orbit 




3 








1 







TABLE II. Types of orbits of massless test particles in the KS space-time. The thick lines represent the range of the orbits. 
The turning points are shown by thick dots. The horizons are indicated by double vertical lines. 



Note that for = we can write the effective potential as 14ff(r) = iV^(r) — 1. This leads to the observation that 
the value of the effective potential at the horizons r± is given by Voff = — 1. Since for black hole solutions we will 
always have Vcsir^)) < — 1 we find that for massive particles r_ < rg < r+. We show the effective potential for £ = 1, 
m = 1 and different values of w in FigO We thus find two different possible radial orbits depending on the value of 
. For — 1 > the particle moving on a radial geodesic will be able to reach the physical singularity at r = 0, 
while that with — l<i?^ — 1<0 cannot reach r — Q and will be deflected at a finite value of r = r-min- Moreover, 



this latter particle cannot reach r = oo and will be deflected at r 



The turning points are at 



min.inax 



with 



' mm, max 



(IV.4) 



Since Vcs{r 
radial m, 
doing so 



mm, max ) 



1 > ^1, we find that T'min < t*- < rp < r_|_ < r-^ 



The particles are thus trapped on 



radial manyworld orbits moving from Tmax to r,„iii and back to Tmax and crossing the horizons in both directions while 



2. Massless test particles 



In Figs. l(d)||l(f)" we show how the effective potential T4ff (r) for a massless test particle (e = 0) changes for different 



values of and uj with m = 1. The potential possesses always two extrema: one maximum, which for uj = oo is 
located at r = 3m and a minimum. The value of this minimum is negative and increases with decreasing lo becoming 
equal to zero in the extremal limit. The existence of a minimum is a new feature as compared to the Schwarzschild 
case. Again, we have an infinite potential barrier at r = 0. Hence, in contrast to the Schwarzschild case we can now 
have three positive, real- valued zeros of f — V^g (r) if E"^ is smaller than the maximum of the potential. The possible 
orbits are a manyworld bound orbit (MBO) on which the particle crosses both horizons with r^m^ < r < Tmax- In 
addition there is an escape orbit (EO) with rmin < r < oo, where the value of rmin fulfills r,nin > r+. These escape 
orbits are very similar to the ones existing in the Schwarzschild space-time. For E'^ larger than the maximum of the 
potential there is only one positive, real- valued zero of £ — Vcs (r) and the particle moves on a two- world escape orbit 
(TEO). The argument that the particle should always cross both horizons for the manyworld bound orbit and the 
two- world escape orbit (TEO), respectively, is similar to the massive case: since Vcf{{r±) = —e = and the turning 
points are given by E'^ — Vcsii") we find that rmin is always smaller than r_ and Tmax is always larger than r+. 
Again, test particles with non-vanishing angular momentum cannot reach the singularity at r = 0. Our results are 
summarized in Table HIl 

The effective potential for radially moving test particles (L^ = 0) is Veff{r) = 0. Hence, all massless test particles 
will reach the singularity at r = on radial geodesies. 



B. Examples of orbits 



In order to find the motion of massive and massless particles in the KS space-time, we have solved the equation 
(|III.6p numerically using the ODE solver of MATLAB that has a 4th order Runge-Kutta method implemented. The 
relative (resp. absolute) errors of the solution are on the order of 10~^^ (10~^^). 



1. Massive test particles 



In Fig|S]we show manyworld bound orbits (MBO) and bound orbits (BO) (region 1 and 2, see Figl3]), respectively, 
for E — 0.99 and Lz = 7.0. In region 2 the test particle moves on a nearly circular orbit with radius much larger 
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FIG. 8. Examples of manyworld bound orbits (MBO) of a massless test particle (e = 0) with E = 0.3, Lz — 1.75 in the 
space-time of a KS black hole with m = 1 and different values of uj. The dashed circles correspond to the two horizons of the 
KS space-time. Note that we are plotting one radial period during which the particle moves twice from rmin to rmax and back 
again. 



than the horizon radii. The shape of the orbit varies only little when changing uj. In region 1, on the other hand, the 
orbit is quite different for w = 10^ as compared to lu = 0.51. For both values of uj, the test particle crosses the two 
horizons in both directions suggesting that these bound orbits are manyworld bound orbits (MBO). Note that this is 
similar to the case of test particles moving in the Reissner- Nordstrom space- time [13, HI] . 

Due to the infinite potential barrier at r = a test particle with non-vanishing angular momentum coming from 
infinity would be reflected at a finite value of r and would not be able to reach r = in the KS space-time. This 
is shown in FiglJl where we give examples of two-world escape orbits (TEO) of a massive test particle with angular 
momentum = A and energy E = 1.8 for different values of uj and m = 1. For all values of w, the particle crosses 
both horizons, but does not reach r = 0, i.e. the particle approaches the KS black hole from an asymptotically flat 
region, crosses both horizons twice and moves away into another asymptotically flat region. 



2. Massless test particles 



As stated above we now have the possibility of manyworld bound orbits (MBO) for massless test particles which 
are not possible in the space-time of a Schwarzschild black hole. In Fig|S]we give examples of manyworld bound orbits 
(MBO) of massless test particles with angular momentum = 1.75 and energy E = 0.3. The qualitative features of 
the orbits are very similar to the massive case. For all values of w the particle crosses both horizons, but due to the 
inflnite potential barrier can never reach the physical singularity at r = 0. Note that bound orbits (BO) of massless 
test particles moving solely outside the black hole do not exist. 

In Figini we give examples of two-world escape orbits (TEO) of massless test particles with angular momentum 
Lz = A and energy E = 1.8. In this case, the test particle encircles the space-time singularity at r = and crosses 
the horizons while doing so. 

In FiglTU]we give an example of an escape orbit (EO) of a massless test particle with angular momentum Lz — 2.2 
and energy E = 0.47 that is deflected by the K S black hole and comes very close to the horizons, but never crosses 
them. This is for m = 1 and uj — 0.51 (see Fig |10(b")| ). For the same values of energy and angular momentum but 
much larger values of uj the test particle would cross the horizons and move on a two-world escape orbit (TEO) (see 
Fig 10(c) ). These orbits should be compared to predictions recently made for massless test particles passing close by 
a Kerr black hole [H. 



C. Observables 



1. Perihelion shift 



The perihelion shift of a bound orbit of a massive test particle in the space-time of a KS black hole can be calculated 
by using pil.4[) . We flnd for the perihelion shift Sip and the period T of the motion of a massive test particle from 
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m= 1,(i} = 0.51,E = 1.80andL =4.00 



m= 1,(i)=1.53,E = 1.8DandL =4.00 



: 1,(i} = 50.00,E = 1.80andL =4.00 
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(c) 



uj = 50.00 



FIG. 9. Examples of two-world escape orbits (TEO) of a massless test particle (e — 0) with E = 1.8, = 4 in the space-time 
of a KS black hole with m = 1 and different values of uj. The dashed circles correspond to the two event horizons of the KS 
space-time. 



L = 2.20, E = 0.47, m = 1 and e = 



1, (0 = 0.51, E = 0.47 and L = 2.20 



0.S1 

1.00x10' 





; 1, (0 = 10000.00, E = 0.47 and L = 2.20 




(a) 



Vefl(r) 



(b) 



EO 



(c) 



TEO 



FIG. 10. Examples of escape orbits (EO) of a massless test particle (e = 0) with E = 0.47, Lz — 2.2 that passes very close by 

For much larger values of lj (here: uj = 10*) the particl e cross es the 
We also show the corresponding effective potential (Fig. 10(a)| . The 



a KS black hole with m = 1 and uu = 0.51 (Fig. |10(b 
horizon on a two-world escape orbit (TEO) (Fig. 10(c 



red dashed line in Fig. 10(a) corresponds to the value of E^, while the dashed circles in Fig. 10(b)pO(c) correspond to the two 
event horizons of the KS space-time. 



''min to r^nax and back again 



Sip ^2 



Lydr 



rV-E2-/(l + L2/r2) 



- 27r , T = 2 



dr 



(IV.5) 



Our results for m = 1 are shown in Fig llli vifhere we give the value of the rate of the perihelion shift Sip/T in 
dependence on uj. In Fig 11(a) we show the perihelion shift for a manyworld bound orbit (MBO), while in Fig 11(b) 



we show that of a bound orbit (BO). We observe that the perihelion shift of the manyworld bound orbit (MBO) is 
much larger than that of the bound orbit (BO). For both types of orbits the perihelion shift increases with increasing 

UJ. 

We can compare the perihelion shift of a bound orbit (BO) in the KS black hole space-time with that in the 
Schwarzschild space-time. Note that the bound orbits of test particles with energy E'^ = 0.9787 and angular mo- 
mentum Lz = 7.00 are nearly circular (see Fig 6(c) and Fig |6(d)| . Hence, it is a good approximation to use the 
perturbative formula for the Schwarzschild space-time which gives the perihelion shift as function of the mass of the 
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central object Ms [27 



(Ms = 67r- 



ms 



GMs 



(IV.6) 



where G is Newton's constant and c is the speed of hght. We have then computed the value of the perihelion shift 
Stp of the bound orbit of a test particle with energy — 0.9787 and angular momentum — 7.00 in the KS 
space-time with m = 1 and several values of w. Setting this value equal to {6ip)s we can find the corresponding mass 
ms that is necessary to obtain the same value of the perihelion shift in the Schwarzschild space-time. We find that 
ms « 1.109 when comparing with the KS space-time for values of w between unity and 10^, i.e. ms does not vary 
much. This leads to the following observation : to have the same perihelion shift in the KS space-time as compared 
to the Schwarzschild space-time we need a smaller mass of the central body. Moreover, the value of rmin (respectively 
Tmiix) is larger (smaller) for a bound orbit (BO) in the KS space-time as compared to a bound orbit (BO) in the 
Schwarzschild space-time with the same value of the perihelion shift. This would be a way to distinguish KS from 
Schwarzschild space-times. This is shown in Fig 12(a)[ where we give the difference Sr^in = '"min,KS ~ ''min.s of the 



minimal radius in the KS space-time rmin^KS and in the Schwarzschild space-time rmin.s as function of uj for two 
different values of m. We also give the value of Sr^ax — ''max.KS ~ ''max.s of the difference of the maximal radius in 
the KS space-time rniax,KS and the maximal radius in the Schwarzschild space-time Tmax.s- Note that the value of 
the radius of the black hole is between m in the extremal limit and 2m in the Schwarzschild limit. For a stellar black 
hole with radius lOfcm, this would correspond to masses between 3.39 solar masses (for w = oo) and 6.78 solar masses 
(for the extremal limit). We observe that the difference decreases with increasing oj (as expected). For increasing m 
both the difference (5r,nin as well as the difference 6r^s.x increase. For this note that Sr^^a.^ is in fact negative and we 
are giving the absolute value here such that the absolute value of Sr^y^^x decreases with increasing m. 



MBO 




,x 10 



BO 




10 



^ « » 



10' 



10 



10 



10' 



(a) 



MBO 



(b) 



BO 



FIG. 11. The value of the perihelion shift per period T as a function of cj for a massive test particle (e = 1) with energy 
= 0.9787 and angular momentum Lz = 7.00 in the space-time of a KS black hole with m = 1. We show the perihelion shift 
per period T for a manyworld bound orbit (MBO) (left) and for a bound orbit (BO) (right), respectively. 



2. Light deflection 

The deflection of light by a KS black hole can be calculated by using ()III.4[) for an escape orbit of a massless test 
particle (e = 0). The light deflection then reads 

oo 

^ = 2 / ^££^=-^, (IV.7) 
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FIG. 12. The absolute value of the difference frmin ~ rmin.KS — '"min,s (5'"max ~ ^max.KS — fmax.s) between the minimal (maximal) 
radius of the bound orbit (BO) of a massive particle in the KS space-time rmin.KS ('"max,Ks) and the minimal (maximal) radius 
of a bound orbit (BO) in the Schwarzschild space-time rmin.s (f'max.s) is shown in dependence on uj for two different values of 
m (left). We also give the difference Srmin for an escape orbit (EO) of a massless test particle (right). In both cases the energy 
of the particle is _E = 0.99 and the angular momentum Lz = 7.00. 



where rmin is the minimal radius of the orbit. Our results for m — 1 are shown in Fig ll31 where we give the value of 
the light deflection in dependence on the impact parameter b, which is equal to Lz/E in the case that the initial value 
of r is equal to infinity. Note that for values of 5(p larger than 2tt the massless test particle first encircles the black 



hole once or several times before going back to infinity. In Fig 13(a) we show the light deflection of the two-world 
escape orbit (TEO), while in Fig |13(b)] we show the light deflection of the escape orbit (EO). 

For the latter case, we observe that the light deflection increases with decreasing impact parameter. This is very 
similar to the Schwarzschild space-time. Lowering the value of the impact parameter further we find that the light 
deflection diverges at a critical value b — 6crit. This critical value depends on uj and decreases with decreasing lu: at 
u! ~ 10'* the value is close to the Schwarzschild value &crit = 5.1961 ~ -\/27, while 5crit = 5.0950 for uj = 2.00 and 
&crit — 4.6937 for lo = 0.50. Lowering the impact parameter even further we find that the light deflection now decreases 
with decreasing impact parameter. This is a new feature as compared to the Schwarzschild space-time, which however 
also exists in the Reissner- Nordstrom case [13, . This phenomenon can be explained when considering the form 
of the effective potential. Lowering the impact parameter b is comparable to flxing Lz and increasing the energy E. 
For small value of E (large values of b) there are three positive real zeros of f — Ves (r) and the corresponding orbits 
are a manyworld bound orbit (MBO) and an escape orbit (EO). Increasing E (decreasing b) we would then find a 
critical value of E for which E'^ is equal to the value of the maximum of the effective potential. This corresponds to 
an unstable circular orbit for which the value of the light deflection diverges. Increasing E (lowering b) further, E^ 
has only one intersection point at positive r with the effective potential and this corresponds to a two-world escape 
orbit (TEO). 

We also flnd that the light deflection for the escape orbit (EO) decreases with decreasing lu, while for the two- world 
escape orbits the dependence on uj depends on the value of the impact parameter. For very small impact parameter, 
the light deflection decreases with decreasing w, while for b close to 6crit it increases with decreasing uj. 

We can again compare with the Schwarzschild case. For large impact parameter b the light deflection in the 
Schwarzschild space-time can be approximated by ^3 

= ^ . (IV.8) 

We have then computed Sip for test particles with — 5.1961 and different impact parameters b = Lz/E in the 
KS space-time and set these values equal to Sip^ to find the corresponding values ms- We find that for w = 2, 
m = 1 we need to choose ms ~ 1.168 for impact parameter 6 = 40 and tos ~ 1.170 for impact parameter b = 29.85, 



15 



30 



Light Deflection 



S-20 



(0 = 0.50 
- (0 = 2.00 
-(0 = 1.00x10' 




(a) 



1.5 2 2.5 3 3.5 4 
impact parameter 

TEO 



4.5 



5.5 



30 



25 



Light Deflection 



S-20 







(0 = 0.50 






(0 = 2.00 


i [ 
' 




(0 = 1.00x10*. 


1 1 
i ' 

\ \ 


L 























(b) 



7 8 
impact parameter 

EO 



10 



FIG. 13. The value of the light deflection of a massive test particle in the space-time of the KS black hole with m = 1 as a 
function of the impact parameter b = Lz/E. We give the light deflection for the two-world escape orbit (TEO) (left) and for 
the escape (EO) (right) for three different values of lj. Note that the light deflection diverges at & = fecrit with bcrit ~ 4.6937 
for uj = 0.50, fecrit = 5.0950 for cj = 2.00 and 6crit = 5.1961 ^ ^27 for a; = lO"*, respectively. 



respectively, to get the same value of the light deflection. For increasing uj the corresponding mg decreases, e.g. 
ms « 1.15 for ui — 10^. The conclusion is very similar to the one in the case of the perihelion shift; to find the 
same value of the light defiection in the KS space-time as compared to the Schvifarzschild space-time the mass of 
the central body has to be smaller. r„iin of the escape orbit (EO) in the KS space-time is larger as compared to an 
escape orbit in the Schwarzschild space-time for the same value of the light defiection. This would be another method 
to distinguish the KS space-time from the Schwarzschild space-time and is shown in Fig |12(b")| where we give the 
difference Sr^in — r^i^^KS — '^inin,s for the escape orbit (EO) of a massless test particle in dependence on uj. Again, 
we observe that 6r^in decreases with increasing uj and increases with increasing m. 



CONCLUSIONS 



In this paper we have studied the motion of massless and massive test particle in the space-time of the KS black 
hole, which is a static, spherically symmetric vacuum solution of Hoi^ava-Lifshitz (HL) gravity. We have taken the 
viewpoint that Hof ava-Lifshitz gravity is essentially a short-distance modification of General Relativity (GR) and have 
used the GR geodesic equation. We observe that there are some new features as compared to the static, spherically 
symmetric vacuum solution of GR, the Schwarzschild solution. For massive test particles we find that next to bound 
orbits there exist manyworld bound orbits on which the test particles cross the two horizons in both directions. For 
massless test particles we can also have manyworld bound orbits, which do not exist in the Schwarzschild case. There 
exist also escape orbits, which are comparable to the ones in the Schwarzschild case as well as two- world escape orbits, 
which are a new feature. Due to an infinite angular momentum barrier, test particles with non-vanishing angular 
momentum can never reach r = - in contrast to the Schwarzschild case where particles that have crossed the event 
horizon unavoidably move to r = 0. Massless test particles moving on radial geodesies will always go to r = 0, while 
massive test particles moving on these geodesies are either trapped on a manyworld radial geodesic if their energy 
E < 1 or they will reach the singularity at r = for E > 1. 

We have also computed the perihelion shift and the light deflection. The rate of the perihelion shift of the manyworld 
bound orbit is much larger than that of the bound orbit and for both orbits this rate decreases with decreasing ui, 
i.e. it is largest for both types of orbits in the Schwarzschild limit. The light deflection increases with decreasing 
impact parameter for escape orbits, but decreases with decreasing impact parameter for two-world escape orbits. 
For escape orbits the light defiection is decreasing for decreasing w, while for two- world escape orbits it decreases 
(resp. increases) for small (large) impact parameter. Approximate methods have been used in several other papers 
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to constrain the value of the parameter wm flTf . Since we believe that constraints from orbits can only be 

obtained for large value of w, i.e. close to the Schwarzschild limit we have not attempted to recompute the constraints 
since we believe that our exact techniques would more or less give the same numbers as those found in 
The aim of this paper has been to solve the geodesic equation exactly and present the complete set of solutions to the 
geodesic equation. 

Recently, the geodesic equation in another Hofava-Lifshitz black hole space-time has been solved analytically in 
terms of hyperelliptic functions fsO] . It seems possible that in some limiting cases of the KS black hole space-time 
considered here, we can also find analytic solutions. This is currently under investigation [26]. 
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